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Abstract: We calculate running coupling corrections for the lowest-order gluon production cross 
section in high energy hadronic and nuclear scattering using the BLM scale-setting prescription. In 
, the final answer for the cross section the three powers of fixed coupling are replaced by seven factors 
of running coupling, five in the numerator and two in the denominator, forming a 'septumvirate' of 
running couplings, analogous to the 'triumvirate' of running couplings found earlier for the small- x 
J> ! BFKL/BK/JIMWLK evolution equations. It is interesting to note that the two running couplings 
in the denominator of the 'septumvirate' run with complex- valued momentum scales, which are 
complex conjugates of each other, such that the production cross section is indeed real. We use our 
lowest-order result to conjecture how running coupling corrections may enter the full fixed-coupling 
fcT-factorization formula for gluon production which includes non-linear small- x evolution. 
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1. Introduction 

With the advent of the LHC era the physics of parton saturation/Color Glass Condensate (CGC) 
[1-28] needs to produce more quantitative predictions with higher accuracy. To achieve this goal 
running coupling corrections to the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner 
(JIMWLK) [13-20] and Balitsky-Kovchegov (BK) [21-25] evolution equations have been calculated 
recently in [29-32] (see also [33,34]). This led to a significant improvement in the comparison of 
CGC predictions with data from heavy ion and deep inelastic scattering (DIS) experiments [35-37]. 
However, to make such comparisons theoretically consistent for observables like hadron production 
in proton-proton {pp), proton-nucleus (pA) and nucleus-nucleus (AA) collisions one needs to also 
include running coupling corrections into the formulas for particle production. 

While an exact analytic formula for gluon production in AA collisions is still not known, we 
do know that in pp and pA collisions gluon production at the level of classical gluon fields and 
leading-ln 1/x nonlinear quantum evolution is given by the fc^-factorization formula [1,38-42]: 



- 1 - 



Here Y is the total rapidity interval of the collision, Cp = (N 2 — 1)/2N C , boldface variables de- 
note two-component transverse plane vectors k = (k 1 ,^), and (f) p , <$a are the unintegrated gluon 
distributions in the proton and the nucleus, respectively, which are defined by [28,38,40] 

4> A (k,y) = I d 2 bd 2 re~ lkr V 2 r N G (r,b,y) (1.2) 



a s (2tt)3 



and 



Mk,y) = [ d 2 bd 2 re- ih - r V 2 n G (r,b,y). (1.3) 

a s (2vrf J 

In Eq. ( |l~2D the quantity N G (r, b, y) denotes the forward scattering amplitude for a gluon dipole of 
transverse size r with its center located at the impact parameter b scattering on a target nucleus 
with total rapidity interval y. N G (r,b,y) can, in general, be found from the JIMWLK evolution 
equation [13-20] . In the large- N c limit it is related to the quark dipole forward scattering amplitude 
on the same nucleus N(r, b, y) by 

N G (r, b, y) = 2 N(r, b, y) - N(r, b, y) 2 , (1.4) 

where N(r, b, y) can be found from the BK evolution equation [21-25] (see also [43] for approximate 
ways of finding N G beyond the large- N c limit without solving the JIMWLK equation). The quantity 
n G (r, b, y) from Eq. ( |1.3| ) is also a gluon dipole amplitude, but taken in a dilute regime, where it is 
found by solving the linear Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation [44,45]. 

In our notation the projectile or target is referred to as a 'proton' if the transverse momenta of 
interest in the problem are much larger than the projectile's (target's) saturation scale, such that 
only linear evolution is needed to describe such a projectile/ target. Conversely, if the saturation 
scale of the projectile/target is comparable to the momentum scales of interest then we refer to this 
projectile/target as a 'nucleus'. With this notation, Eq. ( |1 . 1| ) is not an assumption, but an exact 
answer which resulted in a non-obvious way from a diagram resummation done in [38,41]. 

As mentioned above, the running coupling corrections have been calculated for the BFKL, 
BK, and JIMWLK evolution equations in [29-32] and for the initial conditions to these evolution 
equations in the Appendix of [46] , yielding a prescription of how to find N G and n G in Eqs. ( |1.2j ) 
and ( |1.3| ) including running coupling corrections. However, there are three more factors of strong 



coupling a s in Eqs. (|1.1|), (|1.2|), and ( |1.3| ) for which one needs to set the scale. The first steps in this 
direction were taken in [46], where it was shown that, in accordance with conventional wisdom, in 
order to include running coupling corrections in Eq. ( |1 . 1| ) one has to slightly re-define the observable: 
on top of gluon production, one has to allow for contributions where the would-be produced gluon 
splits into a collinear gluon-gluon (or quark-anti-quark) pair. To ensure that the particles in the 
pair are really collinear, one can put a bound on the virtuality of the gluon splitting into a pair, by 
requiring that k 2 < A^ oll , where /c M is the four-momentum of the would-be produced gluon and A co u 
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is some collinear infrared (IR) cutoff. In [46] it was demonstrated that, for the inclusive production 
cross section of gluons and collinear GG or qq pairs with the invariant mass less than A co u, at least 
one factor of the coupling is a s (A^ oll ). While it may seem tempting to use this result and replace 
a s in Eq. ( |1.1| ) by a s (A^ oll ) constructing a guess for the final answer, it is not even clear a priori 
that Eq. retains its fixed-coupling, A/r-factorized form after the running coupling corrections 

are included. Moreover, there are factors of fixed couplings in Eqs. ( |1.2| ) and ( |1.3| ) for which one 
also has to specify the momentum scale. Therefore the problem of setting the scale of the running 
couplings for the gluon production cross section needs to be addressed by a detailed calculation, 
which we will perform here. 1 

In this paper we will tackle this problem by calculating the running coupling corrections to 
the lowest order (O(af)) gluon production cross section [50-53]. Our lowest-order result will allow 
us to conjecture how Eq. (|1.1|) would be modified due to running coupling corrections. We will 
calculate the running coupling corrections using the scale-setting prescription due to Brodsky, Lep- 
age, and Mackenzie (BLM) [54]. This prescription assumes that the scale of the coupling set by 
the introduction of resummed single-quark-loop corrections is the scale of the coupling set by the 
full calculation including both quark and gluon loop corrections to the propagators and vertices. 
After finding the Nf corrections due to resummed quark loops, one 'completes the beta function' 
by replacing the Nf in the corrections with — 67r/3 2 - (Since we will be resumming powers of a s Nf 
the calculation will not employ the large- N c limit.) 

The paper is structured as follows. We begin in Sec. [| by briefly reviewing the fixed-coupling 
lowest-order calculation of the high energy inclusive gluon production cross section in quark-quark 
scattering. We then include running coupling corrections in Sec. |3|. We begin with the easier 
bremsstrahlung diagrams in Sec. |3.1| , and then move on to the case of diagrams with the triple- 
gluon vertex in Sec. |3.2| . The final result for the gluon production cross section is given by Eqs. 
( |3.31| ) and (|3.32|) in Sec. |3l| We conclude by conjecturing the running coupling generalization of 
Eq. (|1.1| ) given by Eq. (|4.7|) in Sec. [O] and by discussing the implications of our result on the CGC 
predictions for particle multiplicity dN / drj in heavy ion collisions as a function of the centrality of 
the collision in Sec. [4.2| . 

2. Brief Review of the Fixed-Coupling Calculation 

Consider the leading-order contributions to the production of a gluon in a high-energy quark-quark 
scattering in standard Feynman perturbation theory, which we will use throughout the paper. 
We set up the problem in the experimentally relevant momentum regime such that one incoming 
quark has very large "+" momentum = (p^, 0, 0), where we use light-cone coordinates (+, — , _L) 

*A related problem of setting the scale of the running coupling in the formula for the heavy quark pair production 
from [47,48] was studied in [49]. However, in [49] the large masses of the heavy quarks lead to a kinematic regime 
different from the one considered here. 
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throughout the paper with the normalization v ± = (v° ± v 3 )/\/2 for a 4-vector u^, and the other 
incoming quark has very large "— " momentum p% = (0,p^~,0). We take p± and P2 to have the 
same order of magnitude, and pf and p^ are large compared to any other momentum scale in the 
problem. This is the eikonal approximation, and, in particular, p\ and are large compared to 
the transverse momentum of the final state particles. Throughout the paper we will work in the 
light-cone gauge, rj • A — A + = 0, with = (0, 1, 0). In this case the three dominant leading order 
Feynman diagrams are shown in Fig. [[[ which we will refer to as the triple-gluon vertex (A) and 
bremsstrahlung (B and C) diagrams. Note that the two bremsstrahlung diagrams that can be drawn 
with an emitted gluon connecting to the lower quark line are suppressed due to eikonality in this 
A + = light-cone gauge: an easy way to see this is to notice that at high energy the gluon emission 
from the lower quark line should come in to leading order with a factor of 7~ in the quark-gluon 
vertex, which would be multiplied by the gluon polarization component e^, which, in turn, is zero 
in A + = light-cone gauge. 




Figure 1: Diagrams contributing to the lowest-order gluon production in quark-quark scattering at high 
energy in the A + = light-cone gauge. The blob in A denotes a triple-gluon vertex while the outgoing 
gluon line in C is visualized as disconnected for clarity (there are only two gluon lines in C). The initial and 
final quark momenta and fundamental color indices are indicated on diagram A, the adjoint color indices 
are shown in diagrams A and B, and the initial and final quark helicities are indicated on diagram C. 



In a process with three final state on-shell particles, there are five unconstrained momentum 
components in the cross section; we choose to take these components to be the components of 
the transverse momentum transferred from the lower quark line, q, and the components of the 
transverse momentum and the plus momentum of the emitted gluon, k and k + , respectively. 

On-shellness of the final particles in Fig. [I] in particular gives for the top quark line = 
(pi — k + q) 2 ~ —1p\ (k~ — q~), such that 
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to eikonal accuracy (i.e. to leading order in inverse powers of the large momenta, pf and p 2 ). 
Similarly imposing the on-shell condition on the outgoing quark line at the bottom of Fig. [I] yields 

q + = (2.2) 

with the same accuracy. 

The light-cone gauge propagator without loop corrections is 



q A + le q z + te 



g„u - — 

v • 1 



(2.3) 



and the outgoing gluon polarization vector e\ obeys r\ ■ t\ = = and k ■ e\ = for the two 
polarizations A = 1, 2; one may then immediately set 

CA = eI= (°'TT' e )- (2 ' 4) 

Exploiting the eikonality of the process, we can approximate 

Ua' 2 (p2 - q) Y u a2 (p 2 ) ~ 2p$6 ff ' attn , (2.5) 

where <j' 2 and <j 2 are the helicities of the outgoing and incoming lower line quarks, respectively [55]. 
A similar expression holds for the top quark line. The anti-commutation relations of the gamma 
matrices along with the Dirac equation yields the useful relation [56] 

u(p) j.f = 2p ■ eu(p); p^u(p) = 2p ■ eu{p). (2.6) 

Finally, we will also find the following type of cancellation useful 

u(p 2 )$u(p 2 ) = u(p 2 - <?)(- {f 2 ~ + f 2 )u(p 2 ) =0, (2.7) 

by the Dirac equation. 

With the above notation and machinery one may almost immediate write down the leading 
order in eikonality result for the bremsstrahlung diagrams, 

iMo~ 8ig3ptP * 



e x .k(W) i/ ^) j/ J <!ai 6^ 2 , (2.9) 
where repeated indices are summed over. A slightly more involved calculation yields 



lMA * ~ (f-q)W €X ' (k ~ ^ ^ W*^>i ^ ^ (2 ' 10) 
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for the triple-gluon vertex diagram to leading order in eikonality. 

Summing the results for A4a, -Mb, and M. c and using the commutation relations for the color 
matrices one finds that the matrix element for gluon production to leading order in fixed coupling 
a s and eikonality is 

k k 

k? 



i{M A + M B + M 



c 



lg 3 ptp 2 



(k-q)\ 

After summing over final states, averaging over initial states, and including the kinematic factors 
one arrives at the following expression for the gluon production cross section [50-53]: 



da 



2a 3 s C F f d 2 q 



d 2 kr dy 



71^ 



,2^2 



E 



e A ■ (k 



e x • k 



k 



k 2 



Summing over gluon polarizations and opening the brackets in Eq. ( j2.12j ) yields 

da 2«jCi? 1 f d 2 q 



7T Z 



k 2 



q 2 (k 



(2.12) 



(2.13) 



d 2 k T dy 

Our goal now is to set the scales for the three couplings in Eq. (|2.13| ) . Before we do that let us note 
that the lowest-order gluon distribution of a quark is (see e.g. [12,40]) 

a s Cf 1 
7i k 2 



<f>(k,y) 



(2.14) 



With the help of Eq. Q2.14Q we can see that Eq. ( |2.13| ) reduces exactly to Eq. ( |1.1| ) for quark-quark 
scattering. 



3. Running Coupling Corrections 

To include running coupling corrections we will follow the BLM scale-setting procedure [54]. We 
will first resum the contribution of all quark bubble corrections giving powers of Nf, with Nf 
the number of quark flavors and a M the physical coupling at some arbitrary renormalization scale 
yU. We will then complete Nf to the full beta-function by replacing 



in the obtained expression. Here 



N f ^-6ir/3 2 (3.1) 



ft = (3.2) 

12 71 



is the one-loop QCD beta-function. After this, the powers of a M P2 should combine into physical 
running couplings 



a s (Q 2 ) = a » Q2 (3.3) 

1 + a M p 2 In ^2- 
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at various momentum scales Q which would follow from this calculation. Throughout this paper 
we will use the MS" renormalization scheme. 

Below we will begin by including running coupling corrections into the bremsstrahlung diagram 
B from Fig. |T] assuming that only a gluon can be produced in the final state. By doing so we 
will not be able to specify the scale for one of the coupling constants in the diagram, obtaining a 
contribution to the cross section that depends on an arbitrary renormalization scale \i. Following [46] 
we will rectify the problem by redefining the gluon production cross-section to include production of 
collinear gluon-gluon and quark-anti-quark pairs with the invariant mass lower than some collinear 
IR cutoff Ap oll . We will see that this new observable will be completely \i- independent and expressible 
in terms of the running coupling constants. 

3.1 Bremsstrahlung Diagrams 

We begin by considering bremsstrahlung diagrams B and C in Fig. It is easier to see how running 
coupling corrections are organized if we consider the production cross section, which is given by 
the square of the sum of the diagrams in Fig. |I[ Consider the square of the diagram B shown in 
the left panel of Fig. |[ Let us include all the quark loop corrections to it which bring in powers 
of a^Nf. The non- vanishing quark loop corrections to diagram B from Fig. [l] squared are shown 
in the right panel of Fig. 0, where we imply that quark loops need to be resummed to all orders 
on the propagators where we inserted several loops. The corrections are limited to quark bubbles 
on the propagator of the gluon line carrying momentum q in the amplitude and in the complex 
conjugate amplitude. The (massless) quark loop corrections on the outgoing gluon line are zero in 
dimensional regularization, since the produced gluon is on mass shell, k 2 = [57]. Throughout the 
paper we assume that the quarks in the loops are massless: quark mass is not needed to fix the 
scale of the coupling constants. 




Figure 2: Inclusion of quark loop corrections to diagram B from Fig. [T] squared. Vertical dashed line 
denotes the final state cut. 
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The result of the resummation of quark bubble corrections to a gluon propagator is well-known 
(see e.g. [56,57]), and this result is particularly simple for the gluon propagator in the light-cone 
gauge that we are working in. With the form of the light-cone gauge propagator without loop 
corrections, Eq. ( |2.3| ), one can straightforwardly show that, if the contribution from each quark 
loop is written as 



[q 2 g a p - Qa qp] i n(g 2 



(3.4) 



then the resummed gluon propagator is 



(qi + ze) [1 



I% 2 )] D " U 



(?)■ 



(3.5) 



One can see that the tensor structure of the bare propagator, Eq. Q2.3| ), is not modified by the loop 
corrections, Eq. (|3.5|). An explicit calculation of n(g 2 ) to leading order due to the quark loop and 
the Nf part of the gluon propagator counterterm in the MS renormalization scheme yields [57] 



(q 2 + ie) 



where 



6 7T 



/i MS e / . 



DM 



(3.6) 



(3.7) 



We introduce a graphical shorthand for this resummed gluon propagator in Fig. |3|. In the same 
figure we also define for later convenience the cut of this resummed gluon propagator. 



(a) 



■mm, + 



+ 



+ 



(b) 



i i 




+ 



+ 



Figure 3: (a) Graphical shorthand for the resummed gluon propagator, the bare gluon propagator dressed 
by an infinite sum of one loop quark loops, (b) Definition of the cut of the resummed gluon propagator, 
defined here for later convenience. Even though the cuts through gluon propagators on lines with massless 
quark loops are zero in dimensional regularization [57] we keep them explicit here for generality. 



Employing Eq. ( |3.6|) in evaluating the right panel of Fig. §, completing Nf to the full beta 
function using Eq. (|3.1|), and using the definition of the running coupling, Eq. (|3.3j), we see that the 
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fixed couplings of diagram B in Fig. [I] squared are modified by loop corrections to 




(3.8) 



We clearly have a problem, since one of the factors of the coupling is still taken at some 
arbitrary renormalization scale \x. This problem cannot be mitigated simply by the inclusion of 
the corrections due to diagram A, as the above result is all one would obtain in QED for photon 
production (bremsstrahlung) . It seems that the scale of the coupling is not fixed uniquely for 
the diagram in Fig. |2|. As was shown in [46] the problem can be remedied if one redefines the 
observable: on top of the gluon production we should allow for the production of collinear gluon- 
gluon (GG) and quark-anti-quark (qq) pairs. This means that instead of calculating the diagram 
on the right-hand-side of Fig. ^], one has to calculate the graph pictured on the left of Fig. |] As is 
shown in Fig. |3], the cut through the resummed gluon propagator in Fig. |] includes the contribution 
from the square of diagram B in Fig. [j] and also cuts through both gluon propagators with the 
gluon line dressed with quark loops and cuts through the quark loops themselves. Even though 
the contribution from the cuts through gluon propagators with quark loop dressings are zero, as 
mentioned above, we will keep them explicit for the sake of generality (it will turn out that they 
are exactly canceled by a contribution from part of the cuts through the quark loops). Because 
of the redefined observable we are calculating, in addition to the gluon production resulting from 
the cuts through gluon propagators we now have to include the collinear qq pair production from 
cuts through the quark loops. Since we are calculating quark loop corrections to sum up powers of 
a M Nf, we do not need to explicitly include collinear gluon pair production in the calculation: it will 
be included when we complete Nf to the full beta-function. (Note that we are only interested in 
the terms with collinear singularities that contribute to the running of the coupling. Also note that 
- by the BLM prescription - by completing to the full beta-function we are including contributions 
from, e.g., the cut gluon loop, and therefore the running coupling contribution coming from collinear 
gluons.) 

To sum the diagrams represented by the diagram on the left side of Fig. |] we note that the sum 
of the cut diagrams can be written as the imaginary part of a dressed gluon propagator, as shown 
by the diagram on the right side of Fig. [|. The contribution of the dressed propagator is calculated 
using Eq. (|3.b|. Including the adjacent quark-gluon vertices one gets 
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Figure 4: Quark loop corrections to the diagram B in Fig. p] for the calculation of the observable where 
one tags on both the gluons and collinear GG and qq pairs. The left side of the equation indicates the cut 
diagrams to be summed. The right side of the equation explains a simple way of performing the summation 
(see text). The fc-line gluon propagator in the lower panel connects to the upper quark line as indicated by 
quark- gluon vertices at its ends, though this part of the diagram is separated from the rest for illustrative 
purposes. 



where, according to the optical theorem, we took the imaginary part of the "forward amplitude" 
and multiplied it by 2. The factor of i results from the — i needed to convert i M into the amplitude 
M and from i 2 = — 1 coming from the quark-gluon vertices. In the eikonal approximation used 
here the indices of the gluon propagator coupling to the upper quark line carrying a large "+" 
component of momentum are \i — v — +. However, we will keep the indices v unspecified, as 
this will make our discussion more general. In arriving at Eq. ( |3.9| ) we have also completed Nf to 
the full beta function, which allowed us to replace the geometric series by a factor of the running 
coupling constant. Eq. ( |3~9|) is illustrated on the right side of Fig. U it is a short hand way of 
taking the imaginary part of the whole diagram, keeping only the cuts going through the dressed 
gluon propagator with momentum k. 

Eq. ( |3.9|) is still somewhat incomplete, since we need to specify what we mean by the collinear 
qq pairs. To remedy this problem we write using the spectral representation (see [58,59] for a similar 
approach) 



A 2 „ 

coll 



(2vr) 4 J J 2A;+(2tt) 4 J J 2A;+(2 7r) 4 ' 

-oo 

where k 2 = k^ k^ 1 is the virtuality of the gluon line. The last step of Eq. ( pMU| ) completes the 
definition of our observable. By limiting the virtuality of the gluon line between and some small 
momentum scale A^ oll > 0, we keep the diagrams with the gluon line in the final state (top line of 
Fig. |3| (b)) and limit the invariant mass of the produced qq pair (bottom line of Fig. || (b)) to be 
smaller than A^ oll , which ensures the collinearity of the pair. (One can show that a small-virtuality 
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qq pair with both quarks on-shell and massless necessarily has the two quarks collinear with each 
other.) A co u plays the role of the minimum momentum scale to be resolved by the "detector". Below 
we will assume that A co n is much smaller than all the momentum scales involved in the problem, 
but is much larger than Aqcd, such that perturbation theory remains applicable. Resummation 
of higher-order collinear emissions would likely lead to a Sudakov form-factor [60], reducing the 
dependence of the cross section on A co n. (A detailed discussion of the role of A coll can be found 
in [46] , where it is argued that A co n plays the role of the factorization scale for the fragmentation 
functions.) 

Using Eq. ( EHUD in Eq. (p]) yields 



A 2 

il coll 



A 



dk 2 



dk + d 2 k\ 
k+ (2tt)< 



Im 



(-* 



2 p-5/3 



It 



k 2 + ie 



(3.11) 



Since we assume that A co n is the smallest perturbative momentum scale in the problem, we are 
interested only in the terms which do not vanish in the A co n — > limit. Since the numerator 
of the gluon propagator D^ u {k) is completely real and does not have any singularities at k 2 = 
(corresponding to the A co u — » limit), we can move D^k) outside the integral over gluon virtuality 
k 2 and outside the Im part sign. This modifies Eq. ( |3.11|) into 



A 



dk + d 2 k 



A 2 „ 

coll 



k+ (2 7T 



4 D, u (k) I dk 2 Im 



a. 



(-* 



2 p -5/3 



k 2 + ie 



(3.12) 



Extracting the imaginary part in Eq. (|3.12|) we obtain 



A 



dk + d 2 k\ 
k+ (2vr)' 



A coll 

D^(k) J dk 2 



-7T 5{k 2 ) a,(0) + p 



it /3 2 a 2 



(l + ^^lng) +vr2/32«2 



(3.13) 



where for the moment we do not need to specify what we mean by a s (0). Integrating over k 2 in 
Eq. ( EU3D yields 



A w 



X 



dk + d 2 k 
k+ (2tt 



jD, v (k) 



-vra s (0) + — 

P2 



arctan 



*P2a s (A 2 con e-^) 



arctan 



7r/3 2 a s (0) 



(3.14) 



We see now that all the factors of got absorbed into running coupling constants at different 
momentum scales. However, since we are aiming to find the scale of one power of the coupling a M , 
we do not have control over higher powers of the coupling a s which enter ( [3.14 ) . Therefore we need 
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to expand Eq. ( |3.14| ) to the lowest order in a s . This gives 



(3.15) 



with the polarization vector given by Eq. (2^) and * denoting complex conjugation. We see that 
in Eq. ( |3.15| ) we obtained all the standard factors for the outgoing gluon line, along with the running 
coupling constant a s (A 2 oll e _5//3 ) . The kinematic factors and polarizations are the same as for the 
fixed coupling diagram B from Fig. [I] squared. We see that the scale of the remaining coupling 
constant in Eq. ( |3.8| ) has now been fixed and is given by A 2 oll e _5//3 in the MS" renormalization 



scheme. 2 (Note that the terms with a s (0) from Eq. ( |3.14| ) got canceled in arriving at Eq. ( |3.15| ) 



thus relieving us of the need to properly define this quantity.) 

Combining our result in ( |3.15| ) with ( |3.8| ) we see that the three fixed- coupling constants in the 
diagram in Fig. [TIB squared become the following running couplings: 



^« s (ALne- 5/3 )a 2 ( g 2 e- 5 / 3 ). (3.16) 



The above analysis can be repeated for the square of the diagram in Fig. |lp and for the cross 
term between diagrams B and C in Fig. [I], m each case leading to the same answer given by 
Eq. ( |3.16| ). We therefore conclude that Eq. ( |3. gives us the scales of the running couplings for 
the "bremsstrahlung diagrams" B and C in Fig. [I] taken by themselves, without the diagram A. We 
now will analyze the running coupling corrections to diagram A of Fig. [TJ. 

3.2 Triple-Gluon Vertex Diagrams 

Quark loop corrections to the triple-gluon vertex diagram from Fig. [I|A squared are shown in Fig. ^|. 
We now have to "dress" both the q and k — q gluon lines with quark bubbles. The outgoing gluon 
line also needs to be dressed like it was done in Fig. [| for bremsstrahlung diagrams: again the cut 
can go either through the gluon line or through a quark loop. On top of that the triple gluon vertex 
may receive a quark bubble correction as well. Note that in principle the cut can go through the 
quark loop correction to the triple-gluon vertex. However, it is easy to show that that contribution 
does not have a singularity when the gluon virtuality is zero, k 2 = 0. Integration of the contribution 
from such a cut over k 2 from to A 2 oll would thus yield an expression proportional to A 2 oll , which 
vanishes in the A co n — > limit. Since, just like for bremsstrahlung diagrams, we are interested in 
the contributions which do not vanish in the A co n — > limit, the cuts through the quark bubble at 
the triple gluon vertex can therefore be neglected. We see that for the outgoing gluon propagator 



2 Note that for quarks with a nonzero but small mass A 2 oll is replaced by 4 m 2 in the argument of a s (with m the 
quark mass). The quark mass will regulate collinear divergences, making it unnecessary to introduce the IR cutoff 
A 2 
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we have to sum over the same cuts as was done in Sec. O for the bremsstrahlung diagrams, as 
pictured in Fig. |j. Since our derivation of the result Q3.16|) of such a summation was valid for any 
values of the propagator indices /x and u, the conclusion applies to the diagrams in Fig. |5] as well: 
the "dressing" of the outgoing gluon propagator along with the adjacent vertices gives us a factor 



of a s (A 



2 p 

coll e 



-5/3 



k — q 






C.C. 



Figure 5: Quark loop corrections to the square of diagram A from Fig. [l]. The dot denotes the usual 
QCD triple-gluon vertex. 

Concentrating on the top left diagram in Fig. || and using (|3.6|) along with ( |3.1|) we see that 
quark loops lead to the following modification of the factors of coupling in it: 

a 2 

e 



n 2 



1 + /3 2 In 



a. 



(k-q) 2 e- 5 / 3 

7 7^ — 

h MS 



1 2 



Of 



(3.17) 



Indeed, since we have analyzed only one diagram we have not fixed the scales of all the coupling 
constants. However, the expression in ( [3.17 ) is an overall factor, present in each of the diagrams in 
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Fig. |5|. As was noted in [46] the effect of the quark loop corrections to the triple gluon vertex would 
be to generate factors of the type 



1 + ctfi P2 In 



-5/3 



1 + otn 02 In 



A1S 



Q%e 



-5/3 



(3.18) 



A1S 



with 



which, after multiplying Eq. ( p. 171 ), would turn 1/aj; in it into 1/ [a s (Q\ e -5 / 3 ) a s (Q 2 B e~ 5 / 3 ) 
some physical momentum scales Qa and Qb to be determined by an explicit calculation. 

Let us calculate the contribution of a quark-loop vertex correction to the diagram A from Fig. [I]. 
The diagram with the one-loop correction is shown in Fig. || Considering the loop itself (in which 
we include only the quark propagators and quark-gluon vertices along the loop) and adding to it the 
contribution from the loop with the particle number in the loop flowing in the opposite direction 
we find 



1 3 £ 

2 9 f 



abc 



N f 



d d i Tr [r(r + v -it) r/Y (r -/)] 



(2 ir) d {I 2 + i e) [{I - q) 2 + i e] [(I + k - q) 2 + i e] 



(3.19) 



where, in preparation for using dimensional regularization, we have switched to d dimensions in the 
integral. Note that, as before, to simplify the algebra, we assume that all quarks in the loops are 
massless, since this assumption does not affect the scale of the running coupling. 

Before we evaluate (|3.19| ) let us make some simplifications due to the kinematics. The quark- 
gluon vertices at the top (carrying momentum pi) and at the bottom (carrying momentum P2) 
quark lines in Fig. |6| are eikonal: the dominant contribution to the vertex on the top comes with a 
Dirac matrix 7 + , while the vertex at the bottom brings in 7". Using these eikonal vertices along 
with Eqs. (2J.) and (|2.2|) we note that the propagator of the k — q gluon line from Fig. ^ can be 



rewritten as 



jk^w D ^ k - q) ^ 



-i (k-q); 



k+ (k-q) 



(3.20) 



Therefore the v index in Eq. ( |3.19| ) is only transverse, as indicated in Fig. [| The propagator of the 



gluon line carrying momentum q in Fig. ||] can be replaced by 



q 2 q 2 



(3.21) 



making the p-index in Eq. (|3.19| ) equal to +. 

All the simplifications in Eqs. (|2.1| ), (|2T2j ), ( |3.20|) , and (|3.21|) are exactly the same as what 
is usually used in arriving at the fixed-coupling expression for the gluon production cross section 
( |2.12| ). Since we are interested in corrections to the fixed-coupling result, we do not need to keep 
all the factors, since they are the same as for the running-coupling case, too: therefore, to simplify 



the algebra we multiply the expression in Eq. ( |3.19| ) only by (k — q)„ from Eq. ( |3.20| ), g + from 
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Eq. ( |3.21| ), and by eHk) from Eq. ( |2.4| ) for the outgoing gluon line in Fig. [| This modifies Eq. ( |3 . 1 9| ) 



to 



-9 3 f abc N f e x Jk) (k-q); 



d d l Tr [YiV + V {V -^)] 



(2vr) rf (l 2 + ie) [(l-q) 2 + ie] [(I + k - q) 2 + i e] 



(3.22) 



The evaluation of the integral in Eq. (|3.22|) 
is carried out in the Appendix |A| with the answer 
given in Eqs. (TA121) , ( [All ), (^H), and (TATlD . 
In arriving at Eq. (|A12|) we have also added the 
iV/- part of the triple-gluon vertex counterterm 
to remove the infinity in Eq. ( |3.22|) . We note 
that the tensor structure of the one-loop correc- 
tion to the triple-gluon vertex is not limited to 
that of the lowest-order (LO) triple-gluon vertex, 
which, in fact, is well-known in the literature (see, 
e.g., [61,62]). For diagram A in Fig. |l] the con- 
tribution of the triple-gluon vertex is given by 
Eq. ( |A9|) and only leads to terms proportional to 
e\ ■ (k — q) as can be seen from Eq. (|A10 ). The 
full one-loop correction (|A12| ) includes other ten- 
sor structures (see Eq. (|A7|)) leading to the ap- 
pearance of e\ ■ k terms. Of course the ultraviolet 
(UV)-divergent part of the one-loop diagram does 
come in with the tensor structure of the lowest- 
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Figure 6: One quark loop correction to the triple 
gluon vertex in diagram A from Fig. |j. The dis- 
connected arrows denote momentum flow direction, 
while the arrow on the quark loop denotes the par- 
ticle number flow. The diagram with the particle 
number flowing in the opposite direction should also 
be included. 



order triple-gluon vertex (|A9|) . However, in ac- 
cordance with the BLM prescription, we need to collect all a M Nf terms to set the scale of the 
coupling constants: therefore the UV-finite term coming with e\ ■ k in Eq. ( |A12|) is also included. 

Since the tensor structure of the LO triple-gluon vertex is modified, we see that the naive 
expectation for the quark loop correction to simply bring in overall factors like those shown in 
Eq. ( |3.18| ) into the square of the diagram in Fig. | is not correct. Instead we see that the square of 
diagram A in Fig. [I], which at the fixed-coupling order contributed 



(Pkr dy 



2a 3 s C F f d 2 q 



E 



e A • (k 



(k 



(3.23) 



to the gluon production cross section (see Eq. ( |2.12|) above), with the help of Eqs. (|3.17|) and ( |A12j ) 



- 15 - 



gets modified to 



2C f 



d 2 k,Tdy 7r 2 



a. 



(Acoiie 



-5/3 



d 2 q a 2 s (q 2 e~ 5 / 3 ) a 2 ({k - qf e" 5 / 3 ) 



Of 



e x -(k- q) (1 + a M 2 L a ) -e x - k a^(3 2 L b 



(k- q y 



(3.24) 



where L a and are defined in Eq. ( |A13|) and Eq. ( |A14|) , respectively. It is clear that Eq. Q3.24| ) 
can be written in terms of the physical running-coupling constants, eliminating the /i-dependence 
completely. However, we will postpone this last step until the next Subsection, in which we will 
collect all the diagrams. 



k — q 




+ 




Figure 7: Quark loop corrections to the interference terms between diagrams A and B from Fig. |]. The 
notation is the same as in Fig. ||[ 

Before we proceed, let us say a few words about the interference graphs between the triple- 
gluon vertex diagram A and the bremsstrahlung diagrams B and C in Fig. [TJ. The quark loop 
corrections resumming powers of Nf for such interference graphs are illustrated in Fig. [7[ As 
was shown above in the discussion of the diagrams in Fig. ||, cuts through the quark loop correction 
to the triple-gluon vertex do not lead to collinearly-singular contributions (i.e. contributions that 
survive in the A^ oll — > limit) and should be discarded. Therefore the treatment of the quark 
loop corrections for the outgoing gluon propagator should be the same for the interference graphs 
in Fig. ^, as it was for diagram A squared in Fig. [| and also for the diagrams B and C squared 
in Fig. |: in each of the cases we obtain a factor of a s (A^ oll e~ 5 / 3 ). The rest of the quark loop 
corrections is included at the amplitude level, as follows from the above analysis. For instance, the 
right-hand side of each graph in Fig. |7| would bring in a s (q 2 e~ 5//3 ) , while the left-hand side would 
give a contribution proportional to 

a s (qe- 5 / 3 ) a s ({k - qf e^ 3 ) e x ■ (k - q) (1 + a, (3 2 L a ) - e x ■ k a, (3 2 L b 



a, 



[k-q) 



(3.25) 
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as can be inferred from the above discussion and from Eq. (|3.24| ) . 

We are now ready to combine all the diagrams together to determine the scales of all the factors 
of the running coupling. 



3.3 The Result 

Adding the contributions of the "dressed" bremsstrahlung diagrams squared with the help of 
Eq. (|3.16|) to the "dressed" diagram A squared in Eq. (|3.24|) , and including all the interference 
diagrams we finally write the expression for the gluon production cross section with the running 
coupling corrections included: 



da 



2C f 



d 2 kx dy 



71 



(A c 2 oll e- 5/3 ) 



d 2 q a 2 (q 2 e" 5 / 3 ) a 2 ((jfe - qf e" 5 / 3 ) 



,2A2 



a. 



e x -(k-q) (1 + a ljL ^ 2 L a ) - e x ■ k a^^U e A • k ( (fc-q) 2 e 5 / 3 



(3.26) 



where, again, L a and Lb are defined by Eq. ( |A13|) and Eq. (|A14|) , respectively. Equation (|3.26|) can 
be written more compactly as 



da 



2C f 



d?k,T dy 



71 



(AL.e- 5/3 ) 



d 2 q a 2 s (q 2 e~ 5 / 3 ) a 2 ((k - q) 2 e~ 5 / 3 ) 



k — q 
[k^qj-- 



1 + «„ /3 2 In 



Q 

Q 2 e-^ 



2\2 



Qf 



MS 



k 

— ( l + a M /3 2 ln 



2 p -5/3 



Q 2 e 



2 



MS 



(3.27) 



if we define momentum scales Q\ and Q2 by 



ln 



Q 



L _ (k-g) 2 ln^-q 2 ln j- q 2 {k _ q) 2 k 2 (k-q) 2 k 2 [(k - q) 2 + q 2 



MS 



(k — q) 2 — q 2 



[{k-q) 2 -q 2 ^ 



2 [{ k - q f - q 2 



(3.28) 



and 



ln 



Ql , (fc-«) 



In 



fc 5 



MS 



MS 



g 2 (fc-9) 2 [q 2 -(fc-q) 2 -2fc 2 ] (fc-q) 2 
2 [(fc _ q)2 _ q2]3 11 q 2 

| q 2 [(fe _ q) 2 _ g2] + fc 2 [(fe _ g) 2 ± g 2 

2 [(fc - q) 2 - q 2 ] 2 



(3.29) 



We would like to point out that, as expected, the scales Qi and Q2 are independent of fij^, as 
follows from Eqs. (p|) and (^29|) . 
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Recasting Eq. ( 3.27 ) in terms of physical running couplings we get 
da 2C F r ^ 2 



d 2 kx dy 



71 



2 a. (Acoii e" 5 / 3 ) J ^ a 2 (q 2 e^ 3 ) a 2 ((fc - «) 2 e^ 3 ) 



k — q 



k 



[k - qy a s {Q\ e- 5 /3) fc2 ag (Q2 e _ 5 /3) 



(3.30) 



which, together with Eqs. ( |3.28 ) and ( ^j.29| ) can be considered the final answer for the gluon pro- 
duction cross-section with the running coupling corrections included. However, Eq. ( |3.30D is some- 
what unsatisfactory. First, it does not explicitly exhibit the symmetry of the problem under the 
q -H- k — q interchange. Such symmetry follows from the + -H- — interchange symmetry of the 
high-energy scattering of identical hadrons/nuclei. We have broken this up-down symmetry in the 
diagrams by choosing the A + = gauge: however, the physical answer should be indeed indepen- 
dent of the gauge choice, and should be symmetric under q -H- k — q. While Eq. (|3.30|) does posses 
such symmetry, it is not explicitly apparent in the form it is written. Second of all, Eq. ( [3.301 ) does 
not look like the fixed coupling cross-section ( |2.13| ) multiplied by the factors of running coupling, 
which is what one would expect from the procedure of setting the scales of the running coupling 
constants. 

Both of these problems are remedied when, after considerable algebra, Eq. (|3.30|) can be recast 
into the following form: 



da 2C F a s (A 2 oll e" 5 / 3 ) 



d 2 q 



a 2 s (q 2 e-^) a 2 s ((fe - qf e~ 5 / 3 ) 



d 2 k T dy 7T 2 k 2 J q 2 {k-q) 2 a s {Q 2 e~ 5 / 3 ) a s (Q* 2 e" 5 / 3 ) 

with the /^yfg-independent momentum scale Q defined by 

Q 2 I . a'ik-a) 1 1 



(3.31) 



In 



1 q 2 (k - qf 
- m - A 



MS 



MS 



Aq 2 {k-q) 2 [{k-q) 2 -q 



21 6 



k 2 [(k - q) 2 - q 2 ] 



21 3 



X 



[(k - q) 2 } 2 - (q 2 )' (k 2 Y + ((fc - q) 2 - q 2 ) 



,2\2 



2\ 2 



,2\2 



+ 2k 2 



-q 2 (k-q) 2 2 (fc 2 ) 2 + 3 [(k-q) 2 -q 2 Y-3k 2 [(k - q) 2 + q 2 ] 



21 2 



(q 2 ) 3 - [(k - q) 2 } 3 
k-qf 



In 



+ i[(k- q) 2 - q 2 f { k 2 [(fc - q) 2 - q 2 ] k 2 [{k - q) 2 + q 2 } - (q 2 ) 2 - [(* - q) 2 ] 



\2l2 



+ q 2 (k - q) 2 ( k 2 [{k - q) 2 + q 2 } - 2 (k 2 ) 2 - 2 [(fc - q) 2 - q 2 } a ) In 



.212 



[k-q] 



x J 2 q 2 (k - q) 2 + 2 k 2 (k - q) 2 + 2q 2 k 2 - (k 2 ) z - (q 2 )' - [(k - q) 



(3.32) 
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Note that the scale Q 2 is complex- valued! The expression under the square root in Eq. ( 3.32|) is 



non-negative: therefore, to obtain the complex conjugate scale Q* from Eq. ( |3.32|) one only needs to 



change the sign in front of the factor of i in it. The cross-section ( |3.31| ) is, of course, real, as it con- 
tains a complex- valued coupling constant multiplied by its conjugate, a s (Q 2 e~ 5 / 3 ) a s (Q* 2 e~ 5 / 3 ). 

The scale Q 2 is explicitly symmetric under the q ■H- k — q interchange: therefore, the cross 
section ( |3.31| ) is also symmetric under q -h- k — q, as expected from the symmetries of this high 



energy scattering problem. Also Eq. ( |3.31| ) clearly looks like the fixed-coupling cross-section ( |2.13| ) 



with three factors of fixed-coupling replaced by the seven running couplings: we will refer to this 
structure as the septumvirate of couplings. 

Eq. ( |3.31| ), along with Eq. ( |3.32| ), are the main results of this work. They provide us with the 
gluon production cross section in high energy quark-quark scattering with the running coupling 
corrections included. 

Let us point out one interesting feature of our result ( |3.31| ) . It is well-known that the q- integral 



in Eq. ( [3.3 1| ) is dominated by regions where either q w or q w k [1, 2]. Due to the q <H- k — q 



symmetry each of these regions contributes equally, so we can concentrate on only one of them. 
Choosing the q ~ region we see from Eq. ( p.32|) that in the q — > limit 

O 2 k 2 1 

ln-^f- = ln-^- + i (3.33) 

(In fact it is easier to see this by starting from Eqs. ( [3.28D , (|3.29|) , and (|3.30|) .) Neglecting for 
simplicity numerical constants such as e 5//3 and 1/2 we approximate Eq. (|3.31|) by 



k 2 

. 2 \ r j„i 



do _ 4 C F a s (A 2 coll ) r dq_ 2 
d 2 k T dy ~ 7T (k 2 ) 2 J q 2 " s 



(q 2 ) . (3.34) 



The integral over q in Eq. (|3.34|) is cut off in the IR by the saturation scale in the case of hadron- 



hadron, hadron-nucleus, and nucleus-nucleus scattering. Eq. ( ^.34[ ) demonstrates that, loosely 



speaking, the exact Eq. ( ^.31[ ) can be interpreted as having only two powers of ot s on top of ot s (A^ oll ), 
with the two factors of a s running at the smaller of the two scales q 2 and (k — q) 2 . 

Imagine that the quark-quark scattering considered here happens within a larger process of 
hadronic (or nuclear) scattering. After integrating over q 2 with the saturation scale Q 2 as the IR 
cutoff, Eq. (|3.34|) becomes 

da _AC F a s (A 2 coll ) a s (k 2 ) a s (Q 2 ) ^ fc^ 



d 2 k T dy 7T (k 2 ) 2 Ql 

which has a symmetric structure of one coupling running at the 'gluon resolution scale' A^ oll , another 
at the IR cutoff Q 2 characterizing the distribution functions, and the third coupling running at the 
high transverse momentum scale k 2 . 
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4. Discussion 



4.1 Conjecture for the Running Coupling Corrected &; T -Factorization Formula 

Above in Eq. ( |3.31| ) we have calculated running coupling corrections to the lowest-order gluon 



production cross section in high-energy quark-quark scattering. Our result can be easily generalized 
to the case of nucleus-nucleus scattering in the McLerran-Venugopalan (MV) model, in which 
each nucleon in the nuclei is assumed to be made of valence quarks only. The generalization is 
accomplished by multiplying the right-hand side of Eq. ( ^.31[ ) by (N c Ai) (N C A 2 ), with A\ and A 2 
the atomic numbers of the two nuclei. Such a generalization would only be valid for large momenta 
k± = \k\ Qsi,Qs2 where Q s i,Q S 2 are the saturation scales of the two nuclei. In this regime 
non-linear saturation effects are not yet important and can be neglected. 

It would be useful though to generalize our result (|3.31|) (i) beyond the MV model, i.e., to 
include small- x evolution [4,44,45] in it and (ii) inside the saturation region, where the nonlinear 
effects are important [13-25]. Even for the fixed coupling constant result for the gluon production 
cross section including the small- x evolution of [13-25] only exists for the case of proton-nucleus 
scattering, i.e., for the case when = Q s2 3> Q s \ = Q p s [38]; Eqs. ( |1.1| ), (|1.2|), and ( |1.3| ), which 



give the gluon production in proton-nucleus scattering, were shown to be valid in [38] only for 
k± 3> Q p s . Still it would be very useful to generalize Eqs. (|1.1|), (|1.2[), and ( |1.3| ) to include running 
coupling corrections. 

To construct our conjecture for this generalization, let us start working in the framework of 
the MV model. The running-coupling corrections to the Glauber-Mueller formula [63] , the forward 
amplitude of a qq dipole on a nucleus in the MV model, was constructed in the Appendix of [46] 
and is given by Eq. (A8) there. Generalizing this result to the case of a GG dipole scattering on a 
nucleus we write 

"1\ o. 1 



N G (r,b,y = 0) = 1 - exp 



-7T OL s 



^) a s (A 2 ) P T{b)r 2 In. 



A 



(4.1) 



(Strictly speaking the BLM analysis used in [46] to obtain the quark dipole amplitude has to be 
modified for the gluon dipole: however, Eq. ( |4. 1| ) is valid at least in the large- N c limit when it could 
be obtained from Eq. (A8) of [46] with the help of Eq. fll.4|) above.) The dipole amplitude given 
by the Glauber-Mueller formula is rapidity-independent [63] and serves as the initial condition for 
the subsequent BK/JIMWLK evolution: this is denoted by putting y — in Eq. (|4.1| ). Above, 
p is the number density of nucleons in the nucleus and T(b) is the nuclear profile function equal 
to the length of the nuclear medium at the impact parameter b, such that T(b) = 2 \/R 2 — b 2 for 
a spherical nucleus of radius R. A is some scale characterizing the nucleon, which, in general, is 
non-perturbative and is comparable to Aq C d- 

Eq. ( |4.1| ) is an approximation of a more exact expression [46] 

N a (r,b,y = 0) = 1 - exp (-1 pT(b) a GGN (r)\ (4.2) 
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with a GGN (r) the cross-section of the gluon dipole scattering on a nucleon calculated at the two- 
gluon exchange level with the running coupling corrections [46] 

a GGN (r) = [ d*h a 2 ^2 e - 5 /3) ^ _ jl-x _ (4 3) 



212 



[I 

A is the IR cutoff for the Z-integral in Eq. ( |4.3|) . To construct the unintegrated gluon distribution 
at the lowest order we expand Eq. (|4.2j ) to the lowest non-trivial order in a GGN and use the result 
along with Eq. (O) in Eq. (Ol) to obtain 



a: 



(fc 2 e- 5 / 3 ) 1 



a ll <p L A °(k,y = 0) = N C A^— '- —. (4.4) 

In arriving at Eq. ( |4.4| ) we have, for simplicity, assumed that the nucleus is a cylinder of radius R 
with the axis parallel to the beam direction and T(b) = 2 R, such that we could replace pT(b) — > 
N c A/(tt R 2 ) with A the atomic number of the nucleus. 

Eq. ( |4.4| ) demonstrates one important point: it shows that the product <p(k, y) can be rewrit- 
ten in terms of running couplings. It may be possible to define the unintegrated gluon distribution 
<p(k, y) in such a way that it would be expressible in terms of running couplings all by itself. How- 
ever, Eq. (|3.31| ) appears to suggest that such a separation is not necessary, as it contains two 
powers of a s (q 2 ) which are likely to be absorbed into one distribution function and two powers of 
a s ((k — q) 2 ) which are likely to enter the other distribution function. 



Defining 



{k,y) = a s <p{k,y) (4.5) 



we use Eq. ( |4.4|) to rewrite Eq. ( 3.3 1| ) (multiplied by AN 2 for the pA case) in the form of a 
generalization of Eq. (|1 . 1|) : 

da 2C F 1 f 2 -lo, n ,-LO, u a s (A 2 coll e- 5 / 3 ) 

— dq<P p (q,0)4> A (k-q,0) - — — \ 77^-^ . (4.6) 



d 2 k T dy 7t 2 k 2 J yrp v ^ ' a s (Q 2 e" 5 /3) « s (Q* 2 e" 5 / 3 ) ' 



-tLO. 



Again it seems natural that all factorized q-dependence enters one distribution function p (q, 0), 

while all factorized (k — q)-dependence enters the other distribution function </>^°(fc — q, 0). The 
terms depending on momenta q and k — q in a way that can not be factorized into separate q- and 
{k — q)-dependent parts are all included in the scale Q, which generates a "vertex correction". 

Eq. (|4.6p is still an exact result as long as all momenta involved are much larger than the 
saturation scale of the nucleus. However, inspired by this formula we would like to conjecture the 
following running-coupling generalization of Eq. ( |1 . 1|) : 

da 2C F 1 f. 2 - . . — . v v a a (A 2 ol i e" 5/3 ) 



d 2 kr dy it 2 k 



J d\ 4> p (q, y) Mk ~q,Y-y) as{Q ^t]°a s (Q* 2 e-^) (4 ' ?) 
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with the distribution functions defined by 

4>a{*,v) = 7S3 / d 2 bd 2 re- ik - r V 2 r N G (r,b,y) (4.8) 

and 

J p {k,y) = 7S3 / V?nG(r,6,y). (4.9) 

Here Nq should be found from the running-coupling BK/JIMWLK evolution [29,30], while nc 
should be obtained from the running coupling BFKL equation [31]. The scale Q is given by 
Eq. (|3.32|) . We hope that Eq. (|4.7|) is valid in the same regime as the original Eq. ( p..l|) : it should be 
true both inside and outside the nuclear saturation region, with and without the non-linear small-x 
evolution. Of course only exact calculations can prove or disprove the ansatz of Eq. (fOj). 

Just like with the BK/JIMWLK nonlinear evolution equations, it is likely that even if Eq. (|4.7| ) 
is valid, obtaining it may require a subtraction similar to the one done in [29, 30] for the running- 
coupling BK and JIMWLK equations, possibly leading, in the end, to an additive UV-finite cor- 
rection. However, similar to the case of BK/JIMWLK equations, it may be that for a choice of 
subtraction which preserves the linear evolution in the running-coupling part of the answer, the 
additive correction would be small [30,32]. Since Eq. ( |4.7|) is likely to correctly take into account 
all the linear running-coupling BFKL evolution, it may still be a good approximation for the exact 
answer, even if there is an additive correction to it. 

4.2 Multiplicity per Unit Rapidity 

Another interesting question to explore is the effect of our main result ( |3.31|) on the integrated 



gluon multiplicity per unit rapidity dN/dy. The latter is related to the hadronic multiplicity in 
pp, pA and AA collisions, as was utilized in [35,64-66]. Indeed to find dN/dy one has to integrate 
Eq. ( |3.31D over all k: at the same time Eq. ( |3.3ip is only valid for k^^> Q s . However, we will use 



the trick which is valid for the fixed-coupling calculations at least in the quasi-classical limit of the 
MV model [67-69]: we will integrate ( |3.31| ) over k± with the saturation scale Q s providing the IR 
cutoff. The result would be proportional to the correct value of dN/dy up to a constant. We have 
to note that the fact that this procedure gives the correct answer for dN/dy (up to a constant) 
does not, in general, guarantee the same in the running coupling case. In this sense our calculation 
below should be treated as an approximation. 

The integration of the exact Eq. ( [3.3 1|) over k appears to be rather daunting: instead we will 



use the approximate formula ( |3.34j ) and write 

00 k 2 

dN „ A 2 , , „ N f dk 2 f dq 2 „ , 9N , 
^ « 4 C F - a, (Ao„ii) J W J ± *l [tf . (4.10) 
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We assume that both nuclei are identical, each with atomic number A: to generalize ( |3.34| ) to the 
case of AA scattering we multiplied it by A 2 (for simplicity we left out the number of valence quarks 
in each nucleon, which would have given us an extra factor of N 2 ). Here S± = n R 2 is the cross- 
sectional area of the nucleus, which is again assumed to be cylindrical. Performing the integrals in 
Eq. (P^D yields 

™. „ 4CfS± (f) 2 ^ii. (4.11, 

To further evaluate Eq. ( [4.1 1| ) , with the goal of determining the A-dependence of dN/dy, one 
has to define the saturation scale Q s . We are working in the quasi-classical framework of the MV 
model. Therefore, the A-dependence of Q s that we are going to find is necessarily limited to the 
quasi-classical regime, and is going to be strongly modified by small- x evolution with the running 
coupling corrections, as shown in [32,70,71]. Even within the classical approximation there appear 
to be two different ways of defining Q s : 



Defining Q s by requiring that the dipole amplitude ( fill) is of the order of one at \r\ = l/Q s , 
or, more precisely, that the exponent in Eq. ( |4.1|) is equal to —1/4 at \r\ = 1/Q S , yields 

Q 2 S =47ra s (Q 2 )a s (A 2 ) A ln ^ (4 .12) 



where we again replaced pT(b) — > A/S±. Using Eq. ( |4.12|) we eliminate factors of A/S± in 
Eq. ([4. 11|) obtaining 



f.^S « £^j> oc-4-. (4.13) 
dy 4vr 2 ln 2 ^ a 2 (A 2 ) In A 

The proportionality follows from the assumption that for large enough Q s we may neglect 

the differences between A and Aqcd in Eq. Q4.12|) and write Q 2 oc A 1 / 3 . Identifying the 



number of participating nucleons N part with A we see that the multiplicity per participant 
(1/Npart) dN/dy, in this naive model with cylindrical nuclei, is a decreasing function of A. 
It seems that running coupling effects alone are not sufficient to describe RHIC heavy ion 
collision data on hadron multiplicity per participant along the lines of [66] . Other effects used 
in [66], like realistic nuclear profiles leading to two different saturation scales Q™ ax and Q™ n 
for the two nuclei appear to be more important in describing the data (see [72,73]). 

(ii) Another, possibly less justified way of defining Q s is to require that 

^(\k\=Q s ,y) = j^-. (4.14) 
(2tt) 2 

This is motivated by the saturation requirement for the unintegrated gluon distribution in the 
IR, though we have to point out that the distributions entering the fc^-factorization formula 
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( PD given by Eqs. Q and ([Of) in fact go to zero in the k T -»■ limit [39,40]. Eq. ( ^14]) 
gives 

02 = 4vr^(Q2) A, (4.15) 



which, when used in Eq. ( [4.1 1|) yields 



d i~^Sf) a - {Al " )xA - (416) 

We see that in this scenario the particle multiplicity per participant does not change with 
centrality, reinforcing our earlier conclusion that running coupling corrections alone are not 
enough to describe the centrality dependence of the multiplicity per participant observed in 
heavy ion collisions at RHIC. 

Indeed our above conclusions are limited to the quasi-classical regime and are derived under 
the assumption that Eq. (|4.11|) , which gives us the multiplicity of particles with k± > Q s , correctly 
describes the centrality-dependence of the total particle multiplicity. More detailed (possibly nu- 
merical) work is needed to derive the centrality dependence of the conjecture ( |4.7| ). Another way 
of including an A dependence in Eqs. ( |4.13| ) and ( f4.16|) is to identify A co n with Q s (or make it 
proportional to some other power of Q s ). However, it is not clear whether this is justified theo- 
retically, since Q s is clearly not the smallest momentum scale in the problem in the full non-linear 
regime. From the phenomenological side such a substitution would introduce an (extra) decrease of 
(1/N part ) dN/dy with A in both Eqs. Q4.13 ) and ( L161 ), making a successful comparison with RHIC 
data more difficult to achieve. 

4.3 Summary 

To summarize, in this paper we have found the running coupling corrections to the lowest order 
gluon production cross section in hadronic and nuclear high energy scattering. Our exact results 
are presented in Eqs. ( |3.31|) and ( p. 32 ). An approximate simplified version of the exact formula is 



given in Eq. ( |3.35| ). Based on the results ( |3.31| ) and ( |3.32| ) we have conjectured the /c T -factorization 
formula for the gluon production cross section with running coupling corrections included, given in 
Eqs. (f4.7|), ( f4.8| ), and ( |4.9| ). We hope future work will verify our conjecture. 
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A. Loop integral evaluation 

The goal of this Appendix is to evaluate the expression in Eq. (|3.22|) . Following the standard 
procedure for evaluating loop integrals in dimensional regularization [56,57] we introduce Feynman 
parameters x and y and shift the integration variable by defining 



p = p - yq » + x (k- q y. 
After dropping the tilde we rewrite Eq. ( |3.22| ) as 

d d l Tr [7" 7 a Y± 1 P 7 + 1 S ] 



(Al) 



1 l-x 



r = ig 3 r bc N f e^(k)(k-q^ J dx J dy J 



2 n) d [P - (1 _ x - y) ( yq 2 + x ( k - q)2)]3 



with 





x {l a h [-V qs-x(k- q) s ] + L h [y qp - % (k - q)p) + Ip l s [y q a + x (k - q) c 
+ [yq a + x(k-q) a ] [yqp-x(k-q)p] [-y q s - x (k - q) s ]} 



x = 1 — x, y = 1 — y 



(A2) 



(A3) 



and a, (3, and S some internal indices. 

The first three terms in the curly brackets of Eq. (|A2|) are evaluated by substituting 



(A4) 



which makes the Dirac traces trivial. Denoting the contribution of these first three terms by T 3 we 
perform the Wick rotation and integrate over I to obtain 



r 3 =-gf abc k + e x -(k-q: 



a » N f 

2tt 



l-x 



dx (1 + x) I dy I In 



[1-x-y) (yq 2 + x (k - q) 



+ 1 
(A5) 



where we have added the Nf piece of the triple-gluon vertex counterterm in the MS renormalization 
scheme. Integrating (|A5|) over x and y and completing Nf to the full beta-function using (|3.1| ) yields 



r 3 = 2gf abc k + a,f3 2 e A • (k - q) 



1 



12 [{k-qf-q 1 



(I9q 2 -16(fe- q ) 2 ) {{k-qf-q 2 ) 



+3(k-q) 2 (A(k-q) 2 -5q 2 ) In 



(k - qf 



+ 3q 2 (4q 2 -3(k-q) 2 ) In 



MS 



MS 



12 



(A6) 
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We now turn our attention to the last term in the curly brackets of Eq. (^2). We denote the 
contribution of this term by T 4 . Wick rotation, integration over I, and evaluation of the trace of 
Dirac matrices, after somewhat convoluted algebra, gives 

1 l-x 



T 4 = g f abc k+ (k - q)i ^ / ,lv / <hj 



2tt 



1 







X 



y q 2 + x (k — q) 2 

[x 2 (k - q) 2 + y (1 + x) q 2 ] - 2 x (1 - 2 x - 2 y) <f q v L + 2 x (1 - 2 x) <f k v L J. (A7) 

Integrating over x and y in ( |A7[ ), replacing iVy with the help of ( |3.1| ), and bringing e^(k) (k — q)^ 
inside the curly brackets we obtain 

'3(k-q) 2 ( q 2 -( k -q) 2 + q 2 \n^) 

4 [(k - q) 2 - q 2 ] 2 +1 
q 2 -(k-q) 2 + q 2 In ^ 



r 4 = 2gf abc k + a,p 2 U x (k-q) 



+ e A • q (k-q) q 



-e x q (k-q) k 



2[(k-q) 2 -q 2 Y 
[(k - qf - q 2 ] [(k -q) 2 + 3 q 2 ] - q 2 [3 (k - qf + q 2 } In ^-f 



2 i(k-q) 2 -q 



21 3 



(A8) 



In order to more efficiently combine Eqs. ( |X6| ) and ( |A8|) we note that the triple gluon vertex in 
the leading-order diagram A in Fig. [I] contributes 



-2 g f abc k + g" u 
such that, when multiplied by e^(k) {k — q)^, it becomes 



2gf abc k + e x -(k-q). 



Therefore, defining 



we finally write for T = T 3 + T 4 



Tlo = 2gf abc k 



(A9) 
(A10) 
(All) 



T = r LO 02 { e A ■ (k - q) L a - e x ■ k L b 



(A12) 



where we define 

(k - q) 2 In 



MS 



-p. 

MS 



(k — q) 2 — q 2 



q 2 (k-q) 2 k 2 ^ (fe - q) 2 | k 2 [(fe - q) 2 + q 2 



2 [(fe - g)2 _ g 2]2 



(A13) 
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and 



q 2 (k-q) 2 [q 2 -{k-qf-2k 2 ] ^ (fc - qf | q 2 [(k - qf - q 2 ] + k 2 \(k - qf + q 



2 [(fc _ q f - g 2 2 _ g)2 _ g2 

(A14) 

Note that the first term on the right-hand side of Eq. ( |A13| ) is formally similar to what was obtained 
for the kernel of the running-coupling BK/JIMWLK evolution equations (see Eq. (86) in [29]). 
However the similarity does not extend beyond this term. 

References 

L. V. Gribov, E. M. Levin, and M. G. Ryskin, Semihard Processes in QCD, Phys. Rept. 100 (1983) 
1-150. 

J. P. Blaizot and A. H. Mueller, The Early Stage of Ultrarelativistic Heavy Ion Collisions, Nucl. 
Phys. B289 (1987) 847. 

A. H. Mueller and J.-w. Qiu, Gluon recombination and shadowing at small values of x, Nucl. Phys. 
B268 (1986) 427. 

A. H. Mueller, Soft gluons in the infinite momentum wave function and the BFKL pomeron, Nucl. 
Phys. B415 (1994) 373-385. 

A. H. Mueller and B. Patel, Single and double BFKL pomeron exchange and a dipole picture of 



high-energy hard processes, Nucl. Phys. B425 (1994) 471-488, |hep-ph/9403256 | 



A. H. Mueller, Unitarity and the BFKL pomeron, Nucl. Phys. B437 (1995) 107-126, 
|iep-ph/9408245|1 . 



L. D. McLerran and R. Venugopalan, Gluon distribution functions for very large nuclei at small 



transverse momentum, Phys. Rev. D49 (1994) 3352-3355, fhep-ph/9311205 



L. D. McLerran and R. Venugopalan, Computing quark and gluon distribution functions for very 



large nuclei, Phys. Rev. D49 (1994) 2233-2241, |iep-ph/9309289 |. 



L. D. McLerran and R. Venugopalan, Green's functions in the color field of a large nucleus, Phys. 
Rev. D50 (1994) 2225-2233, [|hep-ph/9402335|l . 



Y. V. Kovchegov, Non-abelian Weizsaecker- Williams field and a two- dimensional effective color 
charge density for a very large nucleus, Phys. Rev. D54 (1996) 5463-5469, |hep-ph/9605446|l. 



Y. V. Kovchegov, Quantum structure of the non-abelian Weizsaecker-Williams field for a very large 



nucleus, Phys. Rev. D55 (1997) 5445-5455, |hep-ph/9701229(| . 



-27- 



J. Jalilian-Marian, A. Kovner, L. D. McLerran, and H. Weigert, The intrinsic glue distribution at 



very small x, Phys. Rev. D55 (1997) 5414-5428, jhep-ph/9606337 



J. Jalilian-Marian, A. Kovner, A. Leonidov, and H. Weigert, The BFKL equation from the Wilson 



renormalization group, Nucl. Phys. B504 (1997) 415-431, Piep-ph/9701284 |. 



J. Jalilian-Marian, A. Kovner, A. Leonidov, and H. Weigert, The Wilson renormalization group for 
low x physics: Towards the high density regime, Phys. Rev. D59 (1998) 014014, |frep-ph/97063"77||. 



J. Jalilian-Marian, A. Kovner, and H. Weigert, The Wilson renormalization group for low x physics: 
Gluon evolution at finite parton density, Phys. Rev. D59 (1998) 014015, friep-ph/ 9709432|| . 

J. Jalilian-Marian, A. Kovner, A. Leonidov, and H. Weigert, Unitarization of gluon distribution in 
the doubly logarithmic regime at high density, Phys. Rev. D59 (1999) 034007, [ fiep-ph/9807462| |. 

A. Kovner, J. G. Milhano, and H. Weigert, Relating different approaches to nonlinear QCD evolution 



at finite gluon density, Phys. Rev. D62 (2000) 114005, flhep-ph/0004014 



H. Weigert, Unitarity at small Bjorken x, Nucl. Phys. A703 (2002) 823-860, |hep-ph/0004044| . 

E. Iancu, A. Leonidov, and L. D. McLerran, Nonlinear gluon evolution in the color glass condensate. 



I, Nucl. Phys. A692 (2001) 583-645, [hep-ph/0011241 



E. Ferreiro, E. Iancu, A. Leonidov, and L. McLerran, Nonlinear gluon evolution in the color glass 



condensate. II, Nucl. Phys. A703 (2002) 489-538, flhep-pn/0109115 



Y. V. Kovchegov, Small-x F2 structure function of a nucleus including multiple pomeron exchanges, 



Phys. Rev. D60 (1999) 034008, jhep-ph/9901281 



Y. V. Kovchegov, Unitarization of the BFKL pomeron on a nucleus, Phys. Rev. D61 (2000) 074018, 
|iep-ph/9905214|1 . 



I. Balitsky, Operator expansion for high-energy scattering, Nucl. Phys. B463 (1996) 99-160, 
[)hep-ph/9509348|1 . 



I. Balitsky, Operator expansion for diffractive high-energy scattering, hep-ph/970641 1 



I. Balitsky, Factorization and high-energy effective action, Phys. Rev. D60 (1999) 014020, 
[|hep-ph/9812311|]. 



E. Iancu and R. Venugopalan, The color glass condensate and high energy scattering in QCD, 
|hep-ph/0303204 . 



H. Weigert, Evolution at small x^: The Color Glass Condensate, Prog. Part. Nucl. Phys. 55 (2005) 
461-565, [ |hep-ph/050l087|1 . 



J. Jalilian-Marian and Y. V. Kovchegov, Saturation physics and deuteron gold collisions at RHIC, 



Prog. Part. Nucl. Phys. 56 (2006) 104-231, [|hep-ph/0505052|1 . 



- 28 - 



[29] 

[30] 

[31] 

[32] 

[33] 

[34] 
[35] 

[36] 

[37] 
[38] 
[39] 
[40] 
[41] 
[42] 
[43] 
[44] 



Y. Kovchegov and H. Weigert, Triumvirate of Running Couplings in Small-x Evolution, Nucl. Phys. 



A 784 (2007) 188-226, |piep-ph/0609090 |. 

I. I. Balitsky, Quark Contribution to the Small-x Evolution of Color Dipole, Phys. Rev. D 75 (2007) 



014001, jhep-ph/0609105 



Y. V. Kovchegov and H. Weigert, Quark loop contribution to BFKL evolution: Running coupling and 



leading-N(f) NLO intercept, Nucl. Phys. A789 (2007) 260-284, jhep-ph/061207l l. 

J. L. Albacete and Y. V. Kovchegov, Solving high energy evolution equation including running 



coupling corrections, Phys. Rev. D75 (2007) 125021, [0704.0612 



M. A. Braun, Reggeized gluons with a running coupling constant, Phys. Lett. B348 (1995) 190-195, 
[ |hep-ph/9408261| ]. 



E. Levin, Renormalons at low x, Nucl. Phys. B453 (1995) 303-333, jhep-ph/9412345 ]. 

J. L. Albacete, Particle multiplicities in Lead-Lead collisions at the LHC from non-linear evolution 



with running coupling, Phys. Rev. Lett. 99 (2007) 262301, [arXiv : 0707 . 2545]. 



J. L. Albacete, N. Armesto, J. G. Milhano, and C. A. Salgado, Non-linear QCD meets data: A global 
analysis of lepton- proton scattering with running coupling BK evolution, Phys. Rev. D80 (2009) 



034031, arXiv:0902.1112 



K. Kutak and A. M. Stasto, Unintegrated gluon distribution from modified BK equation, Eur. Phys. 
J. C41 (2005) 343-351, fcep-ph/0408117l] . 



Y. V. Kovchegov and K. Tuchin, Inclusive gluon production in dis at high parton density, Phys. Rev. 



D65 (2002) 074026, [[hep-ph/01 11362 



M. A. Braun, Inclusive jet production on the nucleus in the perturbative QCD with N(c) — > oo, Phys. 
Lett. B483 (2000) 105-114, ^ep-ph/0063003|] . 

D. Kharzeev, Y. V. Kovchegov, and K. Tuchin, Cronin effect and high-p(t) suppression in p a 



collisions, Phys. Rev. D68 (2003) 094013, [|hep-ph/0307037 l. 

Y. V. Kovchegov and A. H. Mueller, Gluon production in current nucleus and nucleon nucleus 



collisions in a quasi- classical approximation, Nucl. Phys. B529 (1998) 451-479, [hep-ph/9802440|. 
A. Kovner and M. Lublinsky, One gluon, two gluon: Multigluon production via high energy evolution, 



JEEP 11 (2006) 083, jhep-ph/0609227 ]. 

Y. V. Kovchegov, J. Kuokkanen, K. Rummukainen, and H. Weigert, Subleading-N c corrections in 



non-linear small-x evolution, Nucl. Phys. A823 (2009) 47-82, [|arXiv : 0812 . 3238 [. 
Y. Y. Balitsky and L. N. Lipatov Sov. J. Nucl. Phys. 28 (1978) 822. 



- 29 - 



E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, The Pomeranchuk singularity in non-Abelian gauge 
theories, Sov. Phys. JETP 45 (1977) 199-204. 

Y. V. Kovchegov and H. Weigert, Collinear Singularities and Running Coupling Corrections to 



Gluon Production in CGC, Nucl. Phys. A807 (2008) 158-189, [|arXiv : 0712 . 3732 



S. Catani, M. Ciafaloni, and F. Hautmann, High-energy factorization and small x heavy flavor 
production, Nucl. Phys. B366 (1991) 135-188. 

J. C. Collins and R. K. Ellis, Heavy quark production in very high-energy hadron collisions, Nucl. 
Phys. B360 (1991) 3-30. 

E. M. Levin, M. G. Ryskin, Y. M. Shabelski, and A. G. Shuvaev, Heavy quark production in parton 
model and in QCD, Sov. J. Nucl. Phys. 54 (1991) 867-871. 

E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, Multi - reggeon processes in the yang-mills theory, 
Sov. Phys. JETP 44 (1976) 443-450. 

Y. V. Kovchegov and D. H. Rischke, Classical gluon radiation in ultrarelativistic nucleus nucleus 



collisions, Phys. Rev. C56 (1997) 1084-1094, [hep-ph/9704201 



A. Kovner, L. D. McLerran, and H. Weigert, Gluon production at high transverse momentum in the 
McLerran-Venugopalan model of nuclear structure functions, Phys. Rev. D52 (1995) 3809-3814, 
Hhep-ph/9505320| j. 

A. Kovner, L. D. McLerran, and H. Weigert, Gluon production from nonAbelian 
Weizsacker-Williams fields in nucleus-nucleus collisions, Phys. Rev. D52 (1995) 6231-6237, 
|hep-ph/9502289(| . 



S. J. Brodsky, G. P. Lepage, and P. B. Mackenzie, On the elimination of scale ambiguities in 
perturbative quantum chromodynamics, Phys. Rev. D28 (1983) 228. 

J. R. Forshaw and D. A. Ross, Quantum chromodynamics and the pomeron. Cambridge University 
Press, 2004. 

M. E. Peskin and D. V. Schroeder, An Introduction to quantum field theory. Addison- Wesley, 
Reading, USA, 1995. 

G. Sterman, An Introduction to quantum field theory. Cambridge University Press, Cambridge, UK, 
1993. 

Y. L. Dokshitzer, D. Diakonov, and S. I. Troian, Hard Processes in Quantum Chromodynamics, 
Phys. Rept. 58 (1980) 269-395. 

Y. L. Dokshitzer and D. V. Shirkov, On exact account of heavy quark thresholds in hard processes, Z. 
Phys. C67 (1995) 449-458. 



- 30 - 



[60] 

[61] 

[62] 

[63] 

[64] 

[65] 

[66] 

[67] 

[68] 

[69] 

[70] 
[71] 

[72] 
[73] 



V. V. Sudakov, Vertex parts at very high-energies in quantum electrodynamics, Sov. Phys. JETP 3 
(1956) 65-71. 

J. S. Ball and T.-W. Chiu, Analytic properties of the vertex function in gauge theories. 2, Phys. Rev. 
D22 (1980) 2550. 

A. I. Davydychev, P. Osland, and L. Saks, Quark mass dependence of the one-loop three-gluon vertex 



in arbitrary dimension, JEEP 08 (2001) 050, jhep-ph/0105072 



A. H. Mueller, Small x Behavior and Parton Saturation: A QCD Model, Nucl. Phys. B335 (1990) 
115. 

D. Kharzeev and M. Nardi, Hadron production in nuclear collisions at RHIC and high density QCD, 
Phys. Lett. B507 (2001) 121-128, [ frucl-th/0012025 l. 



D. Kharzeev and E. Levin, Manifestations of high density QCD in the first RHIC data, Phys. Lett. 
B523 (2001) 79-87, [|tiucl-th/0108006| 1. 



D. Kharzeev, E. Levin, and M. Nardi, The onset of classical QCD dynamics in relativistic heavy ion 



collisions, Phys. Rev. C71 (2005) 054903, [hep-ph/0111315 



T. Lappi, Wilson line correlator in the MV model: relating the glasma to deep inelastic scattering, 



Eur. Phys. J. C55 (2008) 285-292, RarXiv: 0711 .3039 ] . 

J. P. Blaizot, T. Lappi, and Y. Mehtar-Tani, On the gluon spectrum in the glasma, Nucl. Phys. 



A846 (2010) 63-82, [|arXiv : 1005 . 0955 



Y. V. Kovchegov, Classical initial conditions for ultrarelativistic heavy ion collisions, Nucl. Phys. 
A692 (2001) 557-582, [^ep-ph/0011252|[ . 



A. H. Mueller, Nuclear A-dependence near the saturation boundary, tiep-ph/0301109] . 



J. L. Albacete, N. Armesto, J. G. Milhano, C. A. Salgado, and U. A. Wiedemann, Numerical 
analysis of the Balitsky-Kovchegov equation with running coupling: Dependence of the saturation 
scale on nuclear size and rapidity, Phys. Rev. D71 (2005) 014003, flhep-ph/0408216| . 



H.-J. Drescher, A. Dumitru, A. Hayashigaki, and Y. Nara, The eccentricity in heavy-ion collisions 



from color glass condensate initial conditions, Phys. Rev. C74 (2006) 044905, [ nucl-th/0605012 ]. 
A. Kuhlman, U. W. Heinz, and Y. V. Kovchegov, Gluon saturation effects in relativistic U + U 



collisions, Phys. Lett. B638 (2006) 171-177, |nucl-th/0604038 l. 



31 



